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Abstract 

In a previous paper, we have shown that a double phi peak structure appears in 
the dilepton invariant mass spectrum if a first order QCD phase transition occurs in 
ultrarelativistic heavy ion collisions. Furthermore, the transition temperature can be 
determined from the transverse momentum distribution of the low mass phi peak. In this 
work, we extend the study to the case that a smooth crossover occurs in the quark-gluon 
plasma to the hadronic matter transition. We find that the double phi peak structure 
still exists in the dilepton spectrum and thus remains a viable signal for the formation of 
the quark-gluon plasma in ultrarelativistic heavy ion collisions. 
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In a previous paper JTJ, we have proposed that a double phi peak structure in the 
dilepton invariant mass spectrum from ultrarelativistic heavy ion collisions can be used 
to confirm the phase transition from the quark-gluon plasma to the hadronic matter. 
Furthermore, the transverse momentum distribution of the low mass phi peak allows us 
to determine the transition temperature. 

The low mass phi peak results from the decay of phi mesons with reduced in-medium 
mass during the transition. The reduction of the phi meson mass in hot hadronic matter 
is a result of the partial restoration of chiral symmetry B |3|, ^ . In normal effective field 
theory, hadron masses do not decrease at finite temperature and/or density || [J. This 
result is, however, incomplete, because it neglects the modification of the vacuum due to 
medium effects. In QCD, the vacuum can be described in terms of expectation values 
of quark and gluon operators, i.e., condensates. The correlation function for current 
operators defined by quark fields can be calculated as a sum of condensates in the deep 
Euclidean region and can be identified as hadron spectrum function in timelike region. 
Values of the correlation function in these two regions are related to each other through 
the dispersion relation. The behavior of condensates at finite temperature and/or density 
can be calculated with appropriate approximation. Using these condensates on one side 
of the dispersion relation and the hadron spectrum function from an effective theory on 
the other side, it has been shown that hadron masses should decrease in the medium in 
order to satisfy the dispersion relation JF]. In the effective field theory, the vacuum effect 
can be also included via the nucleon antinucleon polarization. Recent studies show that 
this indeed leads to a reduction of meson masses at finite temperature and/or density 
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If a first order phase transition occurs between the quark-gluon plasma and the 
hadronic matter in heavy ion collisions as assumed in Ref. 0, the system spends a 
relatively long time in the mixed phase (about 10-15 fm) during which the phi meson 
mass stays almost constant at a value different from that in free space. Since the du- 
ration of the mixed phase is not negligible compared to the lifetime of the phi meson 
in vacuum (~ 45 fm), a low mass phi peak besides the normal one thus appears in the 
dilepton spectrum. As the transverse flow during the mixed phase is not appreciable, the 
transverse momentum distribution of the low mass phi meson is largely determined by 
the temperature of the mixed phase and provides thus the information on the transition 
temperature. 

In the scenario described in Ref. |]], we have ignored the following effects: (i) the 
collisional broadening of the phi meson width due to its interaction with hadrons; (ii) the 
increase of the phi meson width in the mixed phase due to its interaction with partons 
in the quark-gluon phase; and (iii) the possibility of a smooth transition from the quark- 
gluon plasma to the hadronic matter instead of the mixed phase. If these effects are large, 
then the secondary phi peak proposed in Ref. fl| may not appear. Effects (i) and (ii) 
have recently been studied in Refs. ]I2| and respectively. It has been found that 
the collision of the phi meson with hadrons increases its width to about 10 MeV while its 
interaction with partons adds another few MeV in the width. The resulting phi meson 
width remains small enough to make the secondary phi peak visible. In this paper, we 
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shall study if the double phi peak structure is still present in the dilepton spectrum in 
the case of a smooth transition from the quark-gluon plasma to the hadronic matter. 
According to recent lattice calculations JL4|], the QCD phase transition is a crossover 



but very close to the first order one. There exists a sudden change in the entropy density 



within a temperature interval of less than ~ 10 MeV |15| . To incorporate these features 
of the QCD phase transition, we parametrize the temperature dependence of the entropy 
density s of the hot matter as follows: 
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where Sh{T) = 12aT 3 and s q (T) = 148aT 3 with a = n 2 /90 are the bag model entropy 
density with two flavors in the hadron and the quark phase, respectively; T c is the critical 
temperature; m, n and V are constants. The typical width of the phase transition is 
given by 2T with T ~ 5 MeV according to lattice calculations. For temperatures that 
satisfy \T — T c \ ^> T, the entropy density given by Eq. ([[]) approaches asymptotically to 
Sh{T) and s q (T), respectively, for T below and above T c . We have introduced m and n to 
include the possibility of an asymmetric phase transition. The case m/n = 1 corresponds 
to a symmetric phase transition. For m/n > 1, the transition is asymmetric, and the 
entropy density changes more in the quark phase than in the hadron phase?} 

The pressure P and the energy density e can be straightforwardly evaluated from 
the entropy density. In this paper, we take the chemical potential to be zero, as we are 
interested in the central region of ultrarelativistic heavy ion collisions, where the baryon 
density is expected to almost vanish. Then, the pressure and the energy density are given 
by 

P(T) = [ T s(t)dt, 
Jo 

e(T) = Ts{T)-P{T). (2) 

We note that only one of the three quantities s, P, and e is independent at fixed tempera- 
ture T and chemical potential. In Figs. 1-3, we show s, P, and e by solid lines as functions 
of temperature. The parameters used in evaluating these quantities are T c = 180 MeV, 
T = 5 MeV, and m/n = 1. The dashed lines show the equation of state obtained from 
the standard bag model. 

In the following calculations, we use the same temperature dependence of the phi 
meson mass as in Ref. 10. This is shown in Fig. 4. The decrease of the phi meson 



■''Strictly speaking, if the phase transition is a crossover, the distinction between the hadron phase and 
the quark-gluon phase does not exist. Here we assume, however, that the system is in the quark-gluon 
phase if T > T c and the hadron phase if T < T c . 



3 



mass is mainly due to the presence of a considerable number of strange particles at high 
temperatures. The details can be found in Refs. |], [ITJ]. As in Ref. U, the rho meson 
mass is taken to have the following temperature dependence 



m p {T) _ r — nV« 



m p (T = 0) 



(3) 



and the omega meson mass is independent of temperature. The difference in the tem- 
perature dependence of vector meson masses is due to the difference in their isospin |~7 



We would like to point out that the double phi peak structure in the dilepton spectrum 
depends more on the existence of the shift of the phi meson mass at finite temperature 
than the magnitude of the shift. 

We assume that in ultrarelativistic heavy ion collisions the system has a cylindrical 
symmetry and is in thermal equilibrium. All formulas used in the following calculation 
are given in Ref. [1J. In particular, we include a normalized smearing function of Gaus- 
sian form for the phi meson mass in order to take into account the experimental mass 
resolution: 

F,(M, m +{T)) = * e -(^WT))72^ (4) 
d v 2tt 



where a is a constant and is taken to be 10 MeV |TI 

Assuming boost invariance, we have carried out a hydrodynamical calculation with 
transverse flow for a central collision of 197 Au + 197 Au. We have modified the code of Ref. 



19| to include the smooth quark-gluon plasma to hadronic transition. We have used the 
following values for the parameters [fLl [2(J: the initial proper time, r = 1 fm, and the 
initial radial velocity at the surface of the cylinder, vq = 0; the initial temperature, T = 
250 MeV, the critical temperature, T c = 180 MeV, and the freeze out temperature, 2/ = 
120 MeV; one half of the typical width of the crossover transition, T = 5 MeV and the 
asymmetry factor of the crossover, m/n = 1. Unless specified otherwise, we use these 
default values in the calculation. In addition, we have assumed that the volume fraction 
of the hadron phase is 1 and at T < T c and T > T c , respectively. 

In Fig. 5, we show by the solid line the invariant mass distribution of lepton pairs 
dN/dMdy. We see that the second phi peak between the omega meson and the normal 
phi meson is still visible as in the case of a first order phase transition [|TJ. The low mass 
phi peak is, however, somewhat broadened as the temperature in the present case does 
not stay exactly at the same value during the transition. For comparison, we have also 
plotted in Fig. 5 the result for an asymmetric case by the dashed line. The asymmetric 
factor is chosen to be m/n = 3. In this case, the entropy density drops much more in 
the quark-gluon phase than in the hadron phase. This makes, however, practically no 
difference in the dilepton spectrum as the temperature stays near the critical temperature 
T c = 180 MeV (doted line in Fig. 6) for a relatively long time in both the symmetric 
and the asymmetric case which are shown in Fig. 6 by the solid and the dashed line, 
respectively. To understand this, we take the Bjorken scaling solution without transverse 
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expansion and denote the critical entropy density as s c , which is given by 

_ ms h {T c ) +ns q {T c ) 
m + n 



(5) 



The duration Ar c for which the temperature remains near constant is then approximately 
given by 



Ar a 
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Since s q (T c ) / Sh(T c ) is very large (= 37/3), Ar c is barely affected by any reasonable change 
of the asymmetry factor. 

In Ref. IPJ, we have pointed out that for a first order phase transition the second phi 
peak does not appear if the initial temperature is less than the critical temperature. This 
is not the case for a smooth crossover as shown Fig. 7, in which the initial temperature is 
taken to be 175 MeV and is less than the critical temperature. The reason for this is that 
for the crossover transition the temperature does not stay exactly at T c but around T c , so 
the second phi peak is expected to be visible if T c — T < T is satisfied even if T < T c . 
The result of our numerical calculation supports this. 

If the critical temperature is small, it is possible that the second phi peak may not 
be observed due to the small separation from the normal phi meson peak. In Fig. 8, we 
show by the solid line the result from the calculation with a lower critical temperature, 
T c = 160 MeV. The second peak is still visible. However, if the width parameter T of 
the crossover transition is larger, the second phi peak becomes broader and merges to the 
normal one. This is shown by the dashed line in Fig. 8 for a larger width parameter, 
T = 10 MeV, and the same critical temperature T c = 160 MeV. We note that as in a first 
order transition our results do not change qualitatively for higher initial temperatures 



21], [22[ P3|1 as long as the critical temperature is kept the same 0. 

In Ref. |]J, we have also pointed out that the critical temperature for the QCD 
phase transition can be extracted reasonably accurately from the transverse momentum 
distribution of the low mass phi meson peak in the dilepton spectrum. A similar proposal 
to measure the phase transition temperature with the transverse momentum distribution 
of dileptons from rho mesons was proposed by Seibert The method with the phi 

meson, however, has several advantages over the latter one: (i) The second peak originates 
exclusively from the matter near the critical temperature, so the signal to noise ratio is 
large; (ii) The transverse flow is still small near the phase transition. Therefore, the 
transverse momentum distribution of the second peak reflects more faithfully the critical 
temperature; and (iii) The phi peak is narrower than the rho peak, which makes the 
subtraction of backgrounds easier. We have repeated the same procedure as we did in 
Ref. and have confirmed that the transition temperature can also be extracted from 
the low mass phi meson peak in the case of the crossover transition. In Fig. 9, we show the 
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slope parameter of the dilepton distribution at small transverse momenta as a function 
of the initial temperature. The solid curve is the slope parameter of the low mass peak 
at about 916 MeV. It changes only slightly as the initial temperature increases. On the 
other hand, the slope parameter of the normal phi meson peak at 1019 MeV changes 
significantly as the initial temperature increases due to the resultant development of an 
appreciable transverse flow in the hadron phase. We note especially the sudden increase 
of the slope parameter of the normal phi meson peak in the region, T c — T<T<T c + r, 
where the entropy density changes rapidly as the temperature increases. 

In summary, due to the reduction of the phi meson mass in a hot matter and the 
sudden change of the entropy density at the phase transition, a distinct low mass peak 
besides the normal one appears in the dilepton spectrum even if the quark-gluon plasma to 
hadronic matter transition in heavy ion collisions is not first order. This concludes a series 
of investigations on the possibility of a double phi peak structure in the dilepton spectrum 
from ultrarelativistic heavy ion collisions. With all factors which could potentially weaken 
the double phi peak structure ruled out, the low mass phi peak is thus a credible tool to 
verify the occurrence and to determine the critical temperature of QCD phase transition 
in future ultrarelativistic heavy ion experiments. 
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Foundation under Grant No. A-1110. 



6 



References 

[1] M. Asakawa and C. M. Ko, Phys. Lett. B322, 33 (1994). 

[2] R. D. Pisarski, Phys. Lett. B110, 155 (1982). 

[3] T. Hatsuda and T. Kunihiro, Phys. Rev. Lett. 55, 158 (1985). 

[4] T. Hatsuda, Nucl. Phys. A544, 27c (1992). 

[5] C. S. Song, Phys. Rev. D49, 1556 (1994). 

[6] M. Hermann, B. L. Friman, and W. Norenberg, Z. Phys. A343, 119 (1992). 

[7] M. Asakawa and C. M. Ko, Phys. Rev. C 48, R526 (1993); 
M. Asakawa and C. M. Ko, Nucl. Phys. A560, 399 (1993). 

[8] S. Gao, R. K. Su, and P. K. N. Yu, Phys. Rev. C 49, 40 (1994). 

[9] H.-C. Jean, J. Piekarewicz, and A. G. Williams, Phys. Rev. C49, 1981 (1994). 
[10] H. Shiomi and T. Hatsuda, University of Tsukuba preprint, UTHEP-276. 
[11] C. S. Song, P. W. Xia, and C. M. Ko, Texas A&M University preprint. 
[12] C. M. Ko and D. Seibert, Phys. Rev. C49, 2198 (1994). 
[13] D. Seibert and C. M. Ko, Kent State University preprint, KSUCNR-002-94. 
[14] F. R. Brown, F. P. Butler, H. Chen, N. H. Christ, Z. Dong, W. Shaffer, L. I. Unger, 

and A. Vaccarino, Phys. Rev. Lett. 65, 2491 (1990). 
[15] N. H. Christ, Nucl. Phys. A544, 81c (1992). 
[16] M. Asakawa and C. M. Ko, Nucl. Phys. A in press. 
[17] T. Hatsuda, Y. Koike, and S. H. Lee, Nucl. Phys. B394, 221 (1993). 
[18] PHENIX Collaboration, Conceptual Design Report, January 1993. 
[19] P. Levai and B. Muller, Phys. Rev. Lett. 67, 1519 (1991). 
[20] M. Asakawa, C. M. Ko, and P. Levai, Phys. Rev. Lett. 70, 398 (1993). 
[21] K. Geiger and B. Muller, Nucl. Phys. B369, 600 (1992); K. Geiger, Phys. Rev. D46, 

4965 (1992); 46 4986 (1992). 
[22] E. Shuryak, Phys. Rev. Lett. 68, 3270 (1992); E. Shuryak and L. Xiong, Phys. Rev. 

Lett. 70, 2241 (1993). 
[23] J. Kapusta, L. McLerran, and D. K. Srivastava, Phys. Lett. B 283, 145 (1992). 
[24] D. Seibert, Phys. Rev. Lett. 68, 1476 (1992). 



7 



Figure Captions 



Fig. 1 Entropy density as a function of temperature. The number of flavor is taken 
to be two and T c is 180MeV. Solid and dashed lines correspond, respectively, to 
the crossover transition with T = 5 MeV and m/n— 1, and the first order phase 
transition with the bag model equation-of-state. 

Fig. 2 Same as Fig. 1 for the energy density. 

Fig. 3 Same as Fig. 1 for the pressure. 

Fig. 4 Temperature-dependent phi meson mass in a hot hadronic matter. 

Fig. 5 Dilepton invariant mass spectrum at central rapidity. The solid curve is the result 
from the hydrodynamical calculations with the default parameter set. The dashed 
curve is obtained by changing the asymmetry factor of the phase transition, m/n, 
to 3. 

Fig. 6 Temperature as a function of proper time. The solid line corresponds to the case 
using the default parameters and the dashed line is obtained by changing m/n to 
3. The dotted line is the critical temperature used in the calculation. 

Fig. 7 Dilepton invariant mass spectrum at central rapidity. The default parameters are 
used except that the initial temperature is T = 175 MeV. 

Fig. 8 Dilepton invariant mass spectrum at central rapidity. For the solid line, the 
default parameters are used except that the critical temperature is T c = 160 MeV. 
For the dashed line, in addition, the width parameter of the transition V is changed 
to 10 MeV. 

Fig. 9 The slope parameter of the phi meson transverse momentum distribution as a 
function of the initial temperature. Except the initial temperature, the default 
parameters are used. Solid and dashed curves correspond to the low mass and the 
normal phi meson peak, respectively. The dotted line is the critical temperature 
used in the calculation. 
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